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Aw PURPOSE Ohms Ure 

This thesis is coneerned with the feasibility Ol eae 
Wind measurements sensed by the SEASAT-A Scatterometer 
(SASS) to derive the surface wind stress boundary conditions 
for use by amy oceanmen1 xed lay ce odcum The SASS) wen 
measurements are a unique source of high spatial resolution 
data that may provide information important to OcCeaniiieaed 
layer processes on otherwise unavailable time angie. 
scales. In addition to high spatial resolution, theese. 
orbital con figuration provided temporal resolution \oraew 
order of 36 hours. A comparison of the number Gi iam 
observations and the number of conventionad Jama 
observations available for a sample one-day period indicates 
the potential impact of this data source (Fig. 1). ge 
potential attributes invite investigation intone 
practicality of using SASS wind measurements to Ime 
ocean model boundary conditions and thus ocean mixed layer 
prediction capability. The quality, quantity and resouniiea 
of the SASS data is addressed with respect to it's potential 
impact on ocean mixed layer predictions. Although thew oure 


not a functional scatterometer currently in Or Di Gye 
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have been developed and are being implemented to fill this 
need. 

Synoptic variability in the oceanic mixed layers 
largely caused by the local atmospheric wind stress and 
surface heat fluxes. The accuracy of the wind speed and it's 
phase with respect to the surface heat flux are critical to 
correctly model ocean mixed layer deepening and shoaling. 
Furthermore, the spatial and temporal resolutionmomeeeme 
surface forcing data determines the space and time scales 
which can be resolved with the mixed layer model. The ocean 
mixed layer model used in this study is a vert iealiy 
integrated, second order closure parameterization (ogee 
turbulent kinetic energy budget (Garwood, 1977). The model 
is one~dimensionalt. Therefore, only the heat ft lux =U 
surface and the entrainment heat flux at the bottom of the 
ned: Vayoe Can cause changes in the temperature profile at 
each geographical location. In the model, these individual 
point profiles are not linked horizontally by any dynamumeal 
mechanisms such as Ekman horizontal advection or Verge 
pumping, Or other interpror fprocesses- Thus, this Va 
focusses on the impact of the SASS winds on UtUhewewe 
dimensional vertical mixing alone. 

For the time period covered by this study, themenme, 
other source of wind data over large Ocean areasmenere 
provided by the Fleet Numerical Oceanography Center (FNOC). 


This data set consists of surface marine winds from the FNOC 
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Northern Hemisphere atmospheric analysis amide prediction 
Povem.  laese Surface Winds are derived from atmospheric 
Merecast model OUtpUL and analysis of irregularly scattered 
observations. 

Primer Mose mmer munis —study 15 to understand the 
differences and similarities between the SASS and FNOC wind 
measurements, and the effect that these differences have on 
ocean mixed layer model forecasts. No study has yet been 
Meorkisned comparing tne FNOC and SASS wind fields. 
Mmerefore, to provide a basis for better explaining model 
comparison results, a direct comparison of the wind fields 
Will also be made. The inherent weaknesses in both wind 
fields estimates will be reviewed in light of the results 


Soeaimned from both direct and model output comparisons. 


B. AN OVERVIEW OF THE STUDY 

tamenapeeretLiethe eharacteristics of the initialization 
Mmielad, and the FNOC and the SASS bevadaryaconalulors used in 
mio Study will be described. Also discussed will be the 
aspects of the SEASAT-A mission which affect the data 
quality. 

HiemeeneoretiCalmand Practical aspects of the space and 
time intervals used for comparison of the FNOC and SASS wind 
Bieiids will be presented in Chapter III. The results of the 
direct comparisons between the wind fields is presented in 


Chapter IV. iemcei ero mes time seme ENOC and SASS wind 


SS 


Fields with an ocean mixed layer model is dlsSCiuls ea. 
Gmapter Vv. A summary of the study, a discussione 7- 
results, and recommendations for the future conc®u@esgae 


Chesis in Chapters ya 


II. DESCRIPTION OF INITIAL AND BOUNDARY CONDITIONS 
A. BACKGROUND 

The model a ae three-dimensional field of ocean 
meimoeravure fOr initialization and verification. The 
Surface heat flux and surface wind stress must bde specified 
Or known functions of time. The ocean thermal structure and 
mies surface heat flux will be discussed briefly, then the 


FNOC and SASS wind fields will be discussed and contrasted. 


ee THE MODEL INITIALIZATION AND VERIFICATION FIELDS 

The model requires a data field which describes the 
initial thermal structure of the upper ocean. The data for 
both this field and the forecast verification field were 
Obtained from the North Pacifie Experiment (NORPAX) TRANSPAC 
Ship-of-opportunity XBT analysis (White and Bernstein, 
1979). The geographic location of the gridpoints for these 
meeldas are depicted as "X"'s in Fig. 2 and coincide with tne 
Anomaly Dynamics Study (ADS) region (Anomaly Dynamics Study, 
1978). The ocean temperatures were obtained by objectively 


analyzing all the XBT profiles for each month at depths of 


mee, 40, 00, 380, 100, 120, 150, 200, 300 and 400m. nS: 
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reduced ADS grid (enclosed by a box). 
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the monthly spatial and temporal scales of the initial- 
maclonetrield are 5 degrees of longitude by 2 degrees of 


latitude. 


of THE SURFACE HEAT FLUX BOUNDARY CONDITION 

The FNOC fields used to define heat flux boundary 
Conditions in earlier studies of the ocean mixed layer have 
proven to be inadequate and have been replaced by 
climatological heat flux fields (Jaramillo, 1984) based on a 


monthly climatology of the North Pacifie (Wyrtki, 1965). 


De THE FNOC SURFACE WIND BOUNDARY CONDITION 

The FNOC wind forcing fields have been derived from 
meecnivea FPNOC Six hourly analysis and prediction fields. 
These fields have been interpolated using a two dimensional 
ror on degree polyoma vOe=nourly Mw time intervals and 
Deo jected onto an 11 by 41 grid (Gallacher, 1979). This 
Pera coincides with the ADS region defined in Fig. 2. 

The spatial scale of the FNOC winds is that of the 63 
by 63 polar stereographic grid (true at 60N) used by the 
PNOC Northern Hemisphere atmospheric model. This spatial 
eeale, ss cnough it Varies with latitude, is adequate to 
resolve the atmospheric synoptic Peeereamane 1Sor Cie order 
of 300km at JON. The temporal scale is 6 hours because 


these fields were archived with this periodicity. (This is 


believed sufficient to resolve synoptic variability faa 


stress). 


Ee THE SASS SURFACE WIND BOUNDARY CONDITION 
1. Background 

The SEASAT satellite carried four remote-sensing 
instruments, in addition to tne Scateeromecer. SEASAT was 
launched into a nearly circular polar orbit (salu. 
approximate altitude of 800 km. SASS was enabled from 1800Z 
Juiliy~6 to 02004 80ctober 10)m 1976 emcee designed to cover 
95 percent of the globe every 36 hours. 

The SASS wind measurements were obtained for this 
study from the Jet Propulsion Laboratory (JPL) Pll GUN 
Data System (PODS). The data records were chosen to coincide 
Withee the AD Saeieted = (Figen) Since only wind speed is 
necessary to derive the wind stress at the sea surface, the 
Che wind direction provided by SASS was neglected tae 
the problem of aliased wind directions, which poSsaagim 
complicates usage of the SASS wind measurements, 1s avoided. 
The wind speed values associated with aliased wind vector 
measurements "usually fall within a few percent of each 
other" (Boggs, 1982), so they were averaged to give one wind 
speed value for each data record. The SASS wind speed 
values have a reference height of 19.5m, as do the FNOC 
Winds, and are calculated assuming an atmospheric boundary 


layer with neutral seach 2 


20 


Zee Walt ty Conv noeneStrlCllons 

Many studies have examined the quality of the SASS. 
data. Jones et al. (1982) and Brown et al. (1982) verified 
miat SASs Wind speeds met the specifications of plus or 
minus 2 m/s or 10% of the wind speed. Hawkins and Black 
(1983) made a verification for near gale force winds of 17- 
18 m/s with similar results. Halberstam (1980) found that 
@eriving wind Speeds at a Standard height reduced the 
effects of stability for an unstable atmosphere. However, 
Brown (1983) pointed to the need for an empirical algorithm 
FOr unstable atmospheric layers which better defined 
mesoscale disturbances, such as fronts. Such an algorithm 
mould allow measurements to be adapted to aapolonen mesoscale 
Meera lons; Otherwise, @€rrors would be produced. Boggs 
@roee) has suggested further constraints on SASS data use, 
maeluding cautious use of wind speed measurements exceeding 
B5en/s. 

With few exceptions, these studies indicate tnat 
the SASS wind speed meaSurementS are aS accurate as 
conventional in situ wind speed measurements when and where 
maey Nave been COmpared. However, caution is necessary in 
meeepving all measured values due to the uncertainty of 
derived results dependent on unknown atmospheric stablity. 
The SASS is able to resolve mesoscale disturbances which 


Mave regions of sharply defined gradients, although the 


a 


absolute values derived for these regions are the most 
uncertain. 

Several quality control restrictions were therefore 
imposed by JPL on the SASS data set. Two of these were used 
Mmitially@iowchis "smudy. 


(1) only wind speeds solutions of 50 m/s or lower were 
accepted as possible aliased solutions; 


(2) only wind speed solutions derived from measurements 
with ineidence angles of 70 degrees or lesS Were 
accepted. 


Examination of the entire data set revealed 54 wind 
speed solutions with a value of 50 m/s. These high wind 
speeds are unlikely for the mid-latitude North Pacis items 
August and suggested the need to impose additional quality 
control measures. The following potential errors in SASs 
measurements were identified: 


(1) Wind speeds measured at low ineidence angles, 20 
degrees or less (i.e., the nadir swath), have a high 
probability of error (Boggs, 1982). The reliapagiee 
of the Bragg scattering relationship decreases 
rapidly as the incidence angle of the radiation 
decreases below 20 degrees (i.e., as it becomes more 
perpendicular to the sea surface). 


(2) Wind speeds measured at high incidence angles tend to 
have erroneously large values. Although igual 
screening allowed incidence angles as high as 70 
degrees, angles above 50 degrees frequently yield 
erroneous wind speed values (Halberstam, 1980 and 
Boggs, 1982). 


(3) The accuracy of wind speeds exceeding 25 m/s is 
questionable because of a lack of ground truth data 
to construct an empirical algorithm for tnose “Spee 
(Boggs, 1982). 


Be 


(4) Separate empirical algorithims for calculating wind 
speeds are used depending on the polarization of the 
Signal from the fore and aft antenna. Tous. wind 
speed measurements derived from one horizontally and 
one vertically polarized signal are more likely to be 
in error (Peter Woiceshyn, personnal communication). 
Of the 268,458 records, 37,767 records (14% of the 
data) were derived with mixed polarization. 
(5) The attenuation coefficient allows the reflected 
microwave signal to be corrected for attenuation due 
YOu TacMOspneric mMOLScture . When the attenuation is 
unknown, Wind speed measurements are subject to error 
from atmospheric moisture (Moore, et al., 1982). The 
total amount of atmospheric moisture was sensed by 
the SMMR (also aboard SEASAT) simultaneously with 
SUriae@e Winds sensed by SASS. The moisture sensing 
was largely restricted to the right swath of the SASS 
path. 
The attenuation coefficient for 136,007 (51%) of 
Che data records was unknown for either the fore or aft 
antenna. For 6,204 records, the forward Antenna Attenuation 
Merrection Coefficient (AACC) was known, and the aft AACC 
was unknown (only 18 of these required a correction to be 
made to the sensed radiation strength). There were no 
measurements with the aft AACC known and the forward AACC 
unknown. The AACC for both antennas was known for 102,236 
records and required no correction to the return microwave 
erenal. The remaining 555 records were corrected for 
attenuation, which was known for both the fore and aft 
antennas. Exclusion of all data with an unknown AACC is 
impractical. However, correcting the wind speed for those 
records with known AACC could improve the results. 


Knowledge of these data set properties resulted in 


moet Oulowing quality control criteria being imposed as an 
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improved screen for records which dare most proba bi wee 
Pre. 


(1) Measurements made with incidence angles of less’ than 
eQ degrees were rejected. 


(2) Measurements made with incidence angles exceeding 50 
degrees were rejected. 


(3) Measurements derived from signals with mixed 
polarization were rejected. 


(4) Wind speed measurements exceeding 25 Wns were 
rejected. 


The box in Table 1 encloses the records rema time 
after imposing criteria (1) and (2). Imposing Criterion 
alone eliminates 37,767 records, 92% of which were filtered 
out by the first two criteria. Criterion (3) deeewiimee 
eliminate any wind speed measurements greater than 25 m/s 
which were not eliminated by the first two criteria. 

All wind speed measurements were calculated using a 
single functional form developed by Wentz (1977), which 
relates wind speed to backKSCatltered irrad#@ance Tiaks 
algorithim was not verified by ground truth for Wind Speed. 
exceeding 25 m/s (Boggs, 1982). There were 15 wind speed 
measurements in excess of 25 m/s which were not eliminated 
by criteria (1), (2) or (3). Table 2 summarizecwie 
characteristics of these 15 measurements and any possible 
association they may have with synoptic events visible from 
GOES imagery. Only three of the fifteen records Convamaen 


known AACC, while the majority of the remaining records were 
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Beottcloned in regions which were high in atmospheric 
moisture. The lack of information necessary to make adequate 
eorrections for attenuation of the SASS microwave signal due 
femeeumOospNeric MoOlsture increases the probability of these 
high wind speeds being erroneous. There is no means by 
which they can be definitively proven to be correct. All 
three of the measurements with a known AACC covered regions 
very likely to have had wind shear zones associatated with 
maemus Or Storm activity. These three measurements were not 
reinforced by other coincident measurements with similar 
Peeves. Ine record for the 22nd of July was coincident with 
(within 0.5 degrees latitude or longitude and within 24 
fmomrs) eit sother records (all having a known AACC) with an 
average wind speed of 10.2 m/s. The record for the 30th of 
July was coincident with eight other records (all with a 
known AACC) having an average wind speed of 15.2 m/s. The 
Pecord for the 6th of August was coincident with eight other 
records (3 with a known AACC), all averaging 16.6 m/s. The 
Maex OF COlncident points to justify the 25 m/s values 
demonstrates the need to impose criterion (4), despite it's 
pemewhat ad hoc nature. 

ies Tigmuile Scum rOUreadd lusmonal qualitwacontrol 
criteria eliminated 76,104 (28 %) of the data records. It 
is believed that the screened SASS data set contains the 
highest quality data which can be derived from the SASS data 


record for this time period. 


ei 


Fide RESOLUTION 


A. BACKGROUND 

The desired resolution for any numerical Mog@eue. 
determined by the physical processes that it is designed to 
Simulate. Determining whether or not these desires can be 
met is a fundamental step in analyzing the usefulness of 
available data sets. In this study, spatial scales Of tne 
ocean mixed layer model are already restricted by the 
Cemperature field (White and Bernstein, 1978) wWeegiiee 
initialize and verify the model. Therefore, spatial scales 
of 2 degrees latitude by 5 degrees longitude are determined 
by the minimum requirements to be met by the data sets 
supplying the Wind fercing boundary sconduemen: 

The temporal resolution of an ocean mixed layer model 
is limited by the temporal resolution of the modems 
boundary condition data sets. The FNOC analyses has wind 
speed values every Six hours, which allows CLemMpigmam 
resolution of waves 24 hours long. The data are spline 
interpolated to intervals of one hour to provide Mogae 
stability, though maintaining a 24-hour resolut Jonmam 
geophysical events. Thus the diurnal cycle is the Shogvere 


period phenomenon resolvable with the FNOC data. 
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Mis ore al isco Mine itemM "Or CNe shoo Gata set was 
@xamined first. Then the temporal resolution of the data 
feaaetQ De determined so that a comparison grid reasonable 
Eoeroum PNOC and SASS data sets could be defined. 

This chapter will review the theoretical scales of the 
dynamic and thermodynamic processes important to ocean mixed 
mayerepredicvions. "Next> =a brief summary of the FNOC wind 
meer esoluvion is presented. meinen, a more detawled survey 
Seecene SASS data resolution characteristics is given. 
Panally, the comparison grid used for the remainder of the 


study is defined. 


Eee THEORETICAL SCALES OF MOTION 

Momentum, mass and heat fluxes through the ocean surface 
meemeresponsible for many of the motions of the ocean. 
Atmospheric pressure variations and winds produce momentum 
fluxes across the air-sea interface causing turbulence, 
Currents, and various types of wave motion. 

A lag in the ocean's response to atmospheric 
merertiations is caused by its relatively high inertia. This 
factor further complicates the understanding of air-sea 
[mrmerackions,  sinee the oceanic response is rarely in 
equilibrium with the changing atmospheric conditions. The 
mee important to mixed layer amoee eee ramees from hours, 


associated with oceanic (inertial) motions and atmospheric 


2g 


diurnal fluctuations, to months, asso@iated with atmosieneare 
seasonal changes. 

The relatively high inertia of the ocean afiee@vewen— 
Spatial scales as well. For example, the scale of synoptic 
eddies is proportional to the Rossby radius of deformaweee 
In the atmosphere this scale is of the order of 1000 km and 
can be resolved with a grid spacing of 250 km. However, in 
the ocean, the Rossby radius of deformation is of he jemeer 
of 100 km and therefore oceanic mesoscale eddaeemae 
resolved only with a grid spacing Of 925 Km Cieecc. 

There are four main scales of atmospheric feta 
events which are associated with various oceanic scales of 
motion. The first includes seasonal atmospheric varuageea 
on the spatial scale of ocean basins. 

On this scale the ocean responds to the atmospheric 
zonal wind and meridional changes in the atmo@Spierame 
temperature gradient. These responses are in the form of 
Ekman pumping/suction, horizontal advection of water naar 
by ocean gyres and changes in the seasonal mixed layer 
depth. 

The second scale is the atmospheric synoOpUtCl@ eee ee 
which encompasses spatial variations of the order of 1000 km 
and temporal fluctuations ranging from one to ten days. The 
ecean respomds to foreiinas 2b this scale with changes dniegime 
mixed layer stability and depth, and formation of @¢@ =a 


DPront se 
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The third scale is the atmospheric mesoscale, which is 
associated with Scour rm@emucewpolar lOWS. ancdemurricanes. 
These events take place at temporal scales of the order of 
hours, and spatial scales ranging from 10 to 100 km. The 
ocean's response at this scale ineludes turbulent storm 
mixing and inertia-gravity wave formation. 

The fourth scale is the diurnal scale produced by the 
daily variation in solar radiation at the sea surface. 
Thermodynamic responses to atmospheric forcing warm (cool) 
the ocean's upper layers. This increases (decreases) the 
hydrostatic stability, which dampens (enhances) the effect 
of turbulent wind mixing. The temporal scale is fixed at 24 
hours. The relevent spatial scale is vertical rather than 
femeezonGal;= 1t iS the Obukoy depth scale, which is 
Gemormined from the ratio of Snear production to bouyancy 
mex Jntco the layer. The diurnal cycle 18s emphasized as a 
PemeoratLe atMospheric foreing scale because 1€ primarily 
affects ocean mixed layer thermodynamics. The ocean mixed 
Pmeeerws response to the diurnal scale forcing may be viewed 
as changes in available potential energy. Radiative cooling 
at night causes convective mixing, which lowers the 
available potential energy and deepens the mixed layer. 
Warming of the photic zone during the day raises the mixed 
layer's available potential energy and increases the layer's 


stability, inhibiting turbulent mixing by the wind. 
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oF RESOLUTION CHARACTERISTICS OF THEME NOGSDs Tia 

As noted in Chapter II, the temporal scale of tne FNOC 
wind data set used for this study 1s 6 hours. )eiggeee. 
Sufficient to resolve the ocean mixed layer's response to 
the avmos phere Mauna lame, cle. 

The spatial resolution of the FNOC data is determined 
by the grid used for the atmospheric model. Ine@ @ritdwinem. 
63 by 63 polar stereographic projection of the Northern 
Hemisphere, true at 60N. Because of the projection, the 
spatial resolution changes with latitude, The s@akee 
however, is clearly synoptic, and varies from 306km at GON 
to 360km at 50N. The ocean thermal structure field has a 2 
degree latitude by 5 degree longitude spatial resolution 
(approximately 240km by 600km). Thus, the FNOC wind fields 
have a finer spatial resolution than the ocean thermal 
structure fields in the zonal direction, but” a egg 
resolution in the meridional direction. Therefore 
the FNOC winds are used as a boundary condition, thGuiees 
Spatial resolution is limited by the ocean thermal StUgieeaie- 


fields in the zonal direction and by the wind fields inmiaae 


Merildiona lwidiseecs rom: 


Be RESOLUTION CHARACTERPSTECS OF THE SASS sha rian 
It iS important to establish an initial overview Of@@ne 
characteristics of the SASS data coverage. Reasonable 


choices must be made in determining a spatial grid scale and 
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Meme interval useful both for comparison with the FNOC data 
Mmeerds and fOr specifying mixed layer model boundary 
fonaltions. 
1. Spatial Coverage 
ier Seometinn son weicmoBASAT-orbDitwerhe effects of 
meewearth's rotation, and the variations in instrument 
Meeraving mode each contributed to SASS's characteristic 
Spatial coverage. While the coverage aspects due to orbital 
Meometry and Garth's rotation were predictable, the 
fluctuations in coverage produced by changes is SASS 
operating modes were random. The following sections 
eSueline the results and implications of the effect that each 
of these three main aspects had on spatial coverage. 
a. The Geometry of the SEASAT Orbit 
SEASAT precessed about 26 degrees of longitude 
every orbit or 365 degrees of longitude every Galy. Pes 
Peerace coverage resulted in a criss-crossing pattern of 
Sieertal paths (Fig. 3). This pattern reveals the inherent 
g@ermemdence of spatial resolution on latitude. There is an 
Meoerease in data density due to satellite swath 
intersections, both between ascending and descending tracks 
and also where convergent tracks overlap as the poles are 
approached. The impact of latitudinal dependence on spatial 
@everagpe Variations, though predictable, is a primary 
meomcerm in resolution considerations. Also, the daily 


precession of the SEASAT orbit by 5 degrees longitude beyond 
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meeomprevious Nadir path yielded a negligible variation in 
@ata coverage across the study area. 
ee Ht toOmeronmmrume Harth's Rotation 

As the earth rotates, it's movement alters the 
m@oeppler frequency shift employed by SASS to track the 
geodetic position of its measurements. Because the antennas 
have a fixed frequency range, the width of the swaths 
covered by the forward antenna is different from that 
Govered by the aft antenna. The antenna that was able to 
Meccive the widest Swath depended upon the direction of the 
Satellite. The forward antenna Swath is widest when 
descending, and the aft swath is widest when ascending. 
imeem, cNOws typical values for the variation incurred by 
this effect. A wind speed cannot be determined unless the 
Meemnd area sensed from the fore and aft antennas overlap. 
Hence, the narrowest swath will determine the width of data 
coverage (Fig. 4). This variation will cause a net increase 
in ground data density with increasing latitude because the 
narrowest swath is wider at higher latitudes. The higher 
the latitude, the wider the swath and the fewer number of 
Beeles needed to cover the area. 

Coe lew onos Operational Modes 

Mitte ci eecu ssecmmedmdircetly irom the various 
data gathering modes in which the scatterometer operated. 
Each mode was produced sit merous COmOinations of 


horizontal and vertical polarization with each of the four 
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antenna sticks (See Table 4). This enabled coverage on 
either the left or right sides of the satellite's pag 
six modes or on both sides,of nadir at once in twomnodee: 
When SASS was operating so that only a single side was being 
sensed, the data coverage was twice as dense as when both 
sides were being sensed. The width of side swaths (Table 3) 
didn't change for various modes. The nadir swath had an 
average width of 140 km when both side swaths were being 
sensed and 90 km when only one side swath was being sensed, 
A gap of 110 km separated the side swaths from the nadir 
swath in all modes. 


d. The Summary of SASS Spatial Coverage 
Characteristics 


As an aid to this summary, the variat momma 
Spatial coverage that result from many of the effects 
discussed above are illustrated in Fig. 5a ance The 
spatial bin size contoured is 1 degree by 1 degree. 

The time increment for each "frame" in Bale 
is 101 minutes (the time required for one satellite orbit). 
Only eight of the 14 time frames for this day are shown. 
The remaining 6 contained no data because the orbit path 
fell outside the study region. As expected, the sensor's 
Spatial coverage, though generally predictable, de@G@awingy 
have a consistent data density. 

In Fig. 5b, the result of combining wa: 


individual orbit passes for a single day is seen. Although 
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Figure 5(b). 


the exact spatial coverage cannot be predicted for aes 
day because of the random changes in operational mode; the 
general character is illustrated by this figure. 

The SEASAT orbital pattern created a dependence 
of spatial coverage on latitude. This is seen in DOU Naa 
5a and 5b by noting the latitudes where various paras 
intersect or form gaps. Also, the precession of Cheweean 
Shown in Fig. 5b, provides excellent spatial coverage as 
long as the operational mode designated allows swaths to be 
sensed on both sides of the nadir path at once. The effects 
of the earth's rotation on spatial coverage are negligible 
over the study region. 

Examination of the SASS data set used initia 
study reveals that spatial coverage for time frames of the 
order of one day or less contain unsystematically large 
data-void regions approximately 10 per cent of the time. 
These regions are oriented north-south and usually extend fo 
the bondeae of the study area. The unpredictable character 
of the coverage is directly attributable to the continual 
changes in the SASS operational mode. 

Evidence of the different operational modes 
activated over the course of a Single day is shown in Fig. 
Sa. Frame 6 of Fig. Sa shows an ascending orbit Wiemann 
enters the area in operational mode ior 2 (sensor swaths on 
both sides of the satellite nadir swath). Part way through 


the frame, the operational mode was switched to a right 
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sided swath (mode 4,6 or 8). This is apparent from the lack 
Saeeddta On the left side of the subtrack and from the 
Pup mine@ Or data Gensity On the right side of the subtrack. 
ine Subtrack in eae@h frame is bordered by data gaps which 
make it distinguishable. 

The data coverage for each of the remaining 
frames may be explained in a similar manner. The 
ineonsistencies in coverage imposed by hundreds of 
relatively rapid operational mode changes (of the order of 5 
minutes in duration) every day are the primary cause of 
meoradic data distribution. The SASS data set may be 
considered as a "worst case scenario" in terms of spatial 
coverage, since the SEASAT mission was aborted long before 
long-term coverage in any one operational mode could be 
Maintained. 

From Siri Hehastratlon, 1e becomes clear tnat 
the SASS data contain adequate spatial coverage to be used 
GO define the boundary conditions for Garwood's mixed layer 
model. 

2. Temporal Resolution 
Having determined that there is adequate spatial 
coverage for an ocean mixed layer model, the temporal 
resolution was examined. The best time scale available for 


comparison of the wind fields was sought. 
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a. Revisit Time 

The time between successive samples of wind at 
a Single ground point is defined as "revisit Cime” Fee 
quantity assumes the instrument operated in a double Swath 
mode and is graphed in Fig. 6 based on predictable coverage 
enaracterisStices 

The interrelationship between spatial and 
temporal coverage is seen by comparing Fig. 3 and Fig. 6. 
The intersection of orbit passes seen in Fig. 3 is related 
to the maximum time required for a revisit as shown in Pio 
6. The mean values of revisit time in Fig. 6 may be 
attributed to the variation of swath width with lative 


b. Dependence of Temporal Resolution on Spatial 
Coverage 


Sinee the Sead area cannot be viewed as a 
Single poimt for temporal resolution considerations yu 
interrelationship between spatial coverage and temporal 
resolution was examined. 

The average number of orbits (measured from 
three randomly selected start times) required to cover the 
following percentage of the study area for fixed spatial 
bin dimensions of 2 degrees latitude by 5 degrees longitude 


was calculated. The results are: 


70% coverage requires 15.3 orbits, 
80% coverage requires 17.6 orbits, 
90% coverage requires 20.0 orbits, 
100% coverage requires 41.6 orbits. 
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Varwactcloummemn wonoo tCemooral resolution witn 
latitude. The mean and maximum time required 
for SASS to re-sense a nadir point varies with 
Maer ude:, 
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The advantage of using an interpolation eno 
to fill gridS Witneesseenan 100epee cent coverage rather 
than decrease temporal resolution accordingly Jeu 
demonstrated by these results. Although these percentage 
calculations are rough averages, nearly a doubling of the 
number of orbits (or halving of temporal resolv i Guy eee 
required to increase the percentage of area covered from 90% 
tOmn) © Ops The disadvantages of using an interpomamaen 
routine to allow an increase in temporal resolution stem 
primarily from the probable error of estimating a "best 
guess" value where measured values are desired. The effects 
of the interpolation scheme on the SASS data fields will be 
discussed 01 Waar er sseC Gren 

Since the 2 degree by 5 degree spatial bin 
dimensions are the limits imposed by the ocean temperature 
fields used by Garwood's model, these results are of 
particular interest. The use of these spatial dimenieaiea 
requires a temporal period of 15.3 orbits (slightly Womee 
than 24 hours) for an area coverage of 70%. These time= 
Space dimensions are acceptable for model use and provide 


adequate coverage for interpolation requirementa. 
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ie MeSoULiS AND DISCUSS ION 

The resolution characteristics of the mixed layer model 
miecrtalization and boundary condition fields, as well as the 
mresce Scales Of motion theoretically important to mixed 
layer modeling have been reviewed. 

The initalization field has horizontal spatial scales 
of about 240km meridionally by 600km zonally. These scales 
will resolve ocean temperature fluctuations with wavelengths 
of motion of the order of 1000km in the meridional direction 
pam 25 00KkKm ian the zonal direction. Treen wee, the 
Mencalization field is able to resolve only ocean events 
teen OCCUr Onmethe atmospheric synoptic scale or larger. 
These events include mixed layer depth fluctuations due to 
Ekman pumping, horizontal advection of water masses by ocean 
gyres and features 1000km or more in width. These 
nhOrizontal scales can also resolve the length of ocean 
Mmeeores, although not their width. The dominant oceanic 
horizontal temperature gradient on these scales is 
meridional. Uimcee tora sowee ne Sdarecciom= or the finest 
MOrizontal resolution, enabling the best available 
aoe tLion of GBhe most important ocean moermat Strucegure 
Meaereacteristics. The vertical scales of the initialization 
field enable resolution of fluctuations with wavelengths of 
80m in the upper 120m of the mixed layer and wavelengths of 


100m below 120m. These vertical scales enable resolution of 
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the seasonal thermocline as well as definition Sota. 
Synoptic mixed Wayemedenr a. 

The spatial resolution of the heat flux (Dome 
condition fields are on the atmospheric SynoODpt1 Gaya 
because they have been derived from climatological data 
which is able to define SST gradient with scalesmor mea 
order 600km. Although, the original spatial resolution of 
the fields developed by Wyrtki (1965) were 2 degrees by 2 
degrees, they were published with a resolution of 10 degrees 
by 10 degrees. The coarser resolution was used f oVrwiaeee 
study. The synoptic heat flux variations resolved by these 
fields provided foreing that did not restrict jean 
resolution for the particular study region and period used. 
The climatological heat flux fields used had a time scale of 
one month. Garwood (1977) showed the importance of 
resolving the diurnal cycle with the heat flux fields to 
ocean mixed layer model dynamics. To provide ieaaems 
resolution, a prediction of the diurnal cycle asSSCGe@ietiweg 
With the radiative heat flux was superimposed on the total 
heat flux field uSing a Sinusoidal curve fit which awa 
function of latitude, julian day and hour of the dajaaea@ 
temporal resolution of the heat flux fields 18S aDiewie. 
define the diurnal cycle. This temporal resolution was gage 
possible by application of theoretical considerations of the 


diurnal heat budget. 
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Miem@esecl lt hom Ot wien i md Tlelds is defined by the 
memospheric model used to compute them. The horizontal 
Peacial scale is of the order of 300km. This will allow 
resolution of atmospheric forcing events with a wavelength 
of 1200km. The temporal scale of the FNOC wind data is 6 
mois, allowing resolution of diurnal forcing fluctuations. 

The SASS wind field resolution is defined by instrument 
marpameters and orbital configuration. Tne wer UZOnta 
Spatial scale is 50km, enabling resolution of mesoscale 
events with fluctuation wavelengths as small as 200km. The 
Bowe cal  Presolution Of Che SASS fields is irregular because 
of sporadic changes in operational modes. However, the 
scale is of the order of 24 hours. This scale is able to 
merorve fluctuations in atmospheric forcing events with 
periods of 4 days. 

These results show that the model haS a Spatial 
MmeroolUCiOn limited by the oceanic thermal structure 
fieectalization field. The FNOC and SASS wind fields have 
finer resolution available, with the SASS fields holding the 
highest potential in this aspect. \ddigronaliv., “operation 
by SASS in modes that enable a continuous double swath 
coverage could reduce the temporal resolution of the SASS 
melds tO a iz2 hour scale. Even so, the FNOC wind fields 
have the highest temporal resolution. Therefore, the model 


meeolucion used for comparison is limited on temporal scales 
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by the SASS wind fields and on spatial scales by 


initialization “fiiseuge 
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IV. DIRECT COMPARISON OF DATA SETS 

A. GENERAL 

Using the results of Chapter III, a grid was chosen to 
Sempare the SASS and FNOC wind fields. The first section of 
Chis chapter reviews the description of the comparison grid 
and the method of reducing the SASS and FNOC data sets onto 
it. Next the methods and results of a direct comparison of 
fe FNOC and SASS data fields are presented. (Wie Se 
meuparilsons provide initial insight into the similarities 
moe dafferences of the fields, and they provide a basis to 


better interpret later ocean model comparisons. 


Be mee COMPARTSON GRID 

The fields are resolved on identical 2 degree latitude 
Pees) Gepree longitude grids. The time interval for both 
mirelds is 24 hours. These space and time scales are a 
compromise between the possible temporal resolution in the 
FNOC data and possible spatial resolution in the SASS data 
which allows the data bases to be sufficiently complete for 
matitvial comparison. A temporal dimension of 24 hours was 
ei@eoaen, Over the FNOC dimension of six hours, to avoid large 
meavial gaps in the SASS data. The Spatial dimensions 


were limited in the zonal direction by the resolution of the 


aa 


NORPAX ADS ocean temperature analyses previously used to 
initalize and verify Garwood's model (Gallacher, et al. 
1983) and in the meridional direction by the Northern 
Hemisphere grid employed by FNOC at the time of the creation 
of the data base, interpolated onto even degrees of 
lacivude. 
VV. Boundarwes of tie Compara soumenm nd 
The comparison grid dimensions are CNOseieae 
correspond to a subset of the ADS grid (designated as the 
"reducedes AD Smeeawra which has been used for previous 
studies (Gallacher et al., 1983). There were two reasons for 
choosing the reduced ADS grid over the ADS grid foOreiaa 
and model comparison of the FNOC and SASS fields: 
(1) the reduced grid may minimize horizontal advective 
effects induced in the temperature field by the 


Kuroshio Current extension; 


(2) areas of sparse ocean data density in the ADS 
region are avoided. 


2. Interpolation of SASS Data onto Comparison Grid 
The SASS data set was truncated in time to cover a 
34 day period from July 16, 1978 through August 18.0 
This particular time frame was chosen tO S)Paniiews 
consecutive mid-month analyses of the upper ocean thermal 
structure field as documented by White and Bernstein (1979). 
The data records for each 24 hour period (centered 


on 1200Z) were sorted into 2 degree latitude by 5 degree 


longitude bins. The centers of these bins are depicted by 


D2 


"*¥"s in Fig. 2. The records for each bin were tnen averaged 
to obtain a single wind speed value for each 24-hour period. 
These steps yielded ai0O by 14 grid of data every 24 hours 
for 34 days. For each 24-hour time frame this grid of data 
was then interpolated, using a two-dimensional fifth order 
Mem ynomial, onto a Similar grid, designated by "“O"™'s in 
meee ce. «LNLS was done to coincide with the ADS grid points 
used for later model simulation. iiee Dp, Odmens. of Chis 
interpolation step were fields ona 9 by 13 grid, every 24 
meurs for 34 days. The interpolated grid was truncated to 
Mme reduced ADS grid (Fig. 2), resulting in a 6 by 12 data 
Macrix for™each of the 34 days. 
pe enverpolabion Of PNOC Data onto Comparison Grid 

The six-hourly FNOC fields were linearly averaged 
in time for each 24 hour period, centered at 1200Z. The 
data set was truncatedin time toa 34 day period beginning 
meee 15, 1978, and ending August 18, 1978. The spatial 
eoverage was then truncated from the ADS grid to the reduced 
Poemerid (Fig. 2). 

Saeco ene ds, ebmemresult was a data file 
Pemeaining 34 daily grids, each dimensioned 6 by 12 and 


covering the reduced ADS grid area. 
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Ge COMPARISON OF WIND FIELD CHARACTER GS TICs 


1. The Spatial Means 


The spatial mean was computed for each 24-nour grid ~- 


for both FNOC and SASS data sets. The resulting time series 
of spatial means is presented in Fig. 7a. The time series 
begins with July 15th, Siewe asmla |. 

No consistently strong bias in wind speed values is 
evident in the time series. In general, the corre@acvaon 
between field means is good, indicating no detectable phase 
shift between data sets. Although the time resolution of 
one day is rather coarse, this characteristic 1S Promieeame 
because it ensures that the same model interaction (pew 
the momentumpeand bouyaney Sillux) boundaw,, conditions is 
maintained regardless of whether the FNOC or the SASS wind 
field is used. Although the spatial means show generally 
good agreement, there is a period of 7 days starting on day 
6 and a period of 4 days starting on day 15 for which the 
SASS means are consistently higher than the FNOC aaa 
Conversely, there is a 4 day period starting with day 20 for 
which the FNOC winds are consistently higher. ieee 
seemingly minor variations (of the order of 1 m/s) in the 
means may, however, be significant for the mixed layer 


model. This will be discussed in detail in Chapters 
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Ze opavidl Standard, Deviations 

Fig. 7b shows the standard deviation for each field 
associated with the means in Fig. 7a. AS with the means, no 
consistently strong bias appears in these values. Thus, the 
SASS and FNOC fields agree not only in the average magnitude 
of the wind measurements, but also in overall spatial 
memeaclon, Cmowen mot necessarily in Spatial patterns. 
This agreement between the overall variance of the fields is 
eeetmportant factor in the prediction of surface wind 
bpapoeos, Since the stress iS proportional to tne square of 
the wind speed. 

freee sceuaen Contours 

Having determined that the means and standard 
deviations are in good agreement, a comparison between the 
isotach patterns of each wind field was made. For modeling 
ocean mixed layer responses to mesoscale atmospheric events, 
Dime morizontal variation of the wind stress is at least as 
important as the average and the variance of the wind 
Stress. Fig. 8 presents a comparison of the percent 
coverage of the SASS data field before interpolation, the 
pattern correlation between SASS and FNOC isotach contours 
and the root mean square (RMS) error between the SASS and 
FNOC fields. 

The coverage of the SASS data was always above 55%, 
where 100% denotes no data gaps for that day. The average 


coverage was 86%. As discussed earlier, the coverage has no 
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Figure 7. (a) Comparison of SASS (solid) and FNOC (dotted 


field mean time series. 
(b) Same as (a) for standard deviation. 


56 


100 


90 


80 


70 


60 


Ore  O-4 Oss 90 


PATTERN CORRELATION PERCENT COVERAGE. 
0 | 


Ss) 


(M/S) 
4 


CQ 


ilo a Ols 


Figure 








~— : 
= 10 15 20 25 a0 Sys 


JO le SEU EBesy) 
Pine =semiccmen (a) Permeent=coveraze of SASS data 
before interpolation. Cem at tern correlation 


Serwecenmoteomeand PNOC isorach analyses. (c) RMS 
error between SASS and FNOC isotach analyses. 


oi 


predictable pattern for individual days. The page 
correlation between FNOC and SASS fields shows very high 
correlation between tne fields for @mesmaga =. The pet teen 
correlation plot is Mirrored by thee Msverronep iia 

There is a strong similarity between the graphsmwer 
percentage of the area covered by SASS data and the pattern 
correlation between SASS and FNOC. This similarity suggests 
that the difference in the isotach patterns on those days is 
related to data gaps in the SASS data field before 
interpolation. Where the coverage percentage is high, the 
pattern correlation and RMS error show very good agreement 
between the SASS and FNOC fields. 

It should be stressed that these contours are 
isotachs and not streamlines. The scale of Variat toma. 
wind speed is considerably smaller than for wind direction, 
making a high correlation between isotachs more diffieum ce 
achieve than for streamlines. 

4.) Dae EXt memeswon. Comse laine 

A brief comparison of the data density contouUmemad 
the isotach patterns for the days with the best and worst 
correlation was made. This was done to determine, in a 
qualitative sense, the conditions that contribute VowGne 
quality of the two data sets. Sea surface pressure analjee. 
revealed the synoptic patterns associated with the FNOC 


isotach patterns. GOES satellite imagery was used as an 
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independent data source to verify synoptic events 
Momtuencing both SASS and =FNOC isotach analyses. 

The Appendix stresses the results for all 34 days 
Meedyin this study, Similar to Figs. 9 and 10, which will be 
presented in the following sections. This allows a complete 
Meview of the details of both SASS and FNOC data sets. 

a. The Best Case 

Day 2 (July 17, 1978) was chosen from Fig. 8 as 
the day on which the data fields compared most favorably. 
ferme o presents the SASS data density contours and isotach 
analyses for both SASS and FNOC., 

The SASS data density (Fig. 9a) is consistently 
meen Over the study region. Mimee Only feplons of- Zero 
coverage are two narrow areas extending from the northern 
border of the ADS region through 40ON and centered zonally at 
174E and 199E. The agreement between the SASS and FNOC wind 
speed fields (Fig. 9b and c) is very good and illustrates 
Che potential reliability of the seatterometer for measuring 
meena Speeds over large areas. The GOES satellite imagery for 
Day 2, compared well with both the FNOC and the SASS fields. 
This imagery indicated the presence of an occluded low 
pressure system centered at 49N and 167W, which accounts for 
the isotach pattern associated with the maximum wind speed 
m@eectOn dominating Fig. 9b and ec. ge. One FNOC surface 
pressure chart for that day verifies the presence of the low 


and shows two high pressure systems, one in the northeast 
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180 190 200 210 220 
LONGITUDE (6) 


(a) SASS data density; July 16, 1978 (Day 2). 
CI = 5, alternating solid and dotted with Oso lem 
(b) SASS isotach analysis; same day. CI = 2, 
format as is (a). (c) Same as (b) for NOG dagen 
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morner Of ene study region and the other in the southwest 
eomner. Sieve rectiratreimterest for modeling iis the 
agreement between the magnitude of the wind speed gradients 
for the two fields. 

Alvmougn the agreement is very good, some 
differences are apparent. Tne FNOC isotach analysiS appears 
tO De shifted east three or four degrees longitude relative 
Zcne Shoo analysis. Ine = position error is most. likely due 
to the use of a model to predict the 6-hourly FNOC winds. 
The model relies on a widely scattered distribution of data 
mem lLeS initial conditions, whereas SASS provides direct 
Soservations with negligible position errors. The 
differences between the fields are so slight that neither 
feemeeoe considered in error without corroborating ground 
brut h measurements, which do not exist. 

Cemite WWOToSa ly Gace 

Dayo 1G) Clu by Siew Gymnas chosen from=@=rieg. 8 
momecic daywon which the data fields compared least 
Ramorably. Pie. SOs Ullicmocmomass big. 9, except for 
Daye 16. 

Mneweo AS Steatameaensity contours showm in 
leeeeroa tl lustrate tne pattern that is characteristic of 
days with many data gaps. iMmewdaga density 1S not 
consistent and it has sharp maxima and minima. The data 
mapo ober ar ee and €xtend vo the study reszion borders. The 


dominant gradients on the atmospheric synoptic scale are in 
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mnie meridional direction. Sinee mesoscale atmospheric 
events also contain strong zonal gradients, the orientation 
of these data gaps iS an important factor in the capability 
to fill gaps accurately by interpolation. This factor will 
be discussed further in the following section. Both the 
GOES imagery and FNOC surface pressure analysis revealed a 
developing synoptic low pressure system whose position 
@oincided with the maximum in the SASS isotach pattern in 
the western half of the study region (Fig. i0b). The FNOC 
isotach analysis showed a wind speed maximum near about 170W 
which corresponded with an increase in the pressure gradient 
between a high pressure ridge extending north along an axis 
as approximately 160W (south of the study region) and the 
developing low. Since there is sufficient data coverage by 
Boo west of 175W and SASS makes direct observations that 
have been shown to be otherwise reliable, the maximum wind 
Speed region portrayed by the SASS isotach pattern is 
believed to be more reliable than the FNOC model analysis. 
However, east of 175W the SASS field has a gap coincident 
with the area of increased pressure gradient. 

The atmospheric Systems associated with the 
isotach maxima on Day 16 are on a smaller scale than those 
pm ay 25 However, these resuits illustrate the advantage 
of using SASS for sensing atmospheric events at and below 


Che synoptic scale. Mier e ye Ore sASS in turn 
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contributes to weak agreement with the lower resolution FNOC 
wind field. 
Ce ~DiSeessionmet Resulvs 

By studying the best and worst cases as well as 
other periods shown in the Appendix, two factors stand out. 
First, reduced pattern correlations and high RMS errors 
primarily result from low percentage of coverage by the SASS 
wind fields before interpolation. The large data gap ama 
worst case SASS field was the primary reason 1 COmiparea 
poorly with the FNOC field, whereas the best case had little 
Or no data gap. The second factor which app Gai@oueme 
contribute to variations in correlation between (them ieiae 
1s the difference between the Spatial resolution @o7 a 
faseddis: The SASS field has a higher spatial res@aiaiwen 
which causes it to show PQAEURRS unresolved by the FNOC 
fields. When the study region is dominated by a SynOpGae 
pattern, the isotach fields c@mpare very Wwele The 
comparison is not good when there are strong mesoscale 
smaller scale disturbances within the study region. The 
ability of SASS to directly measure winds lends credibility 
to the accuracy of the higher resolution SASS winds when and 
where they are sensed. 

In addition, examination of the fHilgimee 
available in the Appendix shows that the orientation of the 
SASS data gaps in relation to the wind speed gradient 


appears to be a significant consideration. Intéerpolaviing 
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across an elongated gap will be more effective when the 
length of any single contour passing through the gap is 
Minimized. This ensures that both curvature changes in the 
eontour and gradient change across the adjacent contours are 
more likely to be resolved. Alternately, interpolation 
across an elongated gap which is ermented parallel to the 
local gradient in the field is more reliable than 
interpolation with the gap oriented perpendicular to the 
local gradient. This is important because the data gaps 
Meescmit in the daily SASS fields are primarily caused by 
mitoses ti Operational mode,;=etnerefore, these gaps are 
eonsistently oriented north-south. Mean wind speed 
gradients on the synoptic scale are generally oriented 
north-south, while mesoscale disturbances may include strong 
east-west gradients. So, when mesoscale disturbances are 
sensed, their orentation to data gaps may cause considerable 
error in the interpolated values. Thus the orientation of 
data gaps and the percentage coverage of the SASS fields 
both degrade SASS's capability to resolve mesoseale events. 
Been Of these factors should influence selection of an 
Setimal number of satellites and their orbit trajectories in 


future observing systems. 


25 


5. Wind S0pessmGemparison 
The wind stress on the ocean is the average 
turbulent transfer of horizontal momentum across thew 


ocean interface by the vertical component of the turbulence, 


T=, (u' w') = jo uxe - ‘al 


Where u' and w' are the variations from the mean 
horizontal and vertical surface velocities, and KP is 
Che density of water. The surface frictionmveloci’, (ane 
the scale of the turbulent velocity in the upper part of the 
ocean mixed layer. The surface stress can be related to the 


mean wind speed using the bulk formula, 


C = (Pal P) Cp UF 2) 
where Cp is the drag coefficient at the height of (item 
Wind speed, u AI oes is the dens iG ae maine. The Saigarw 
coefficient at the sea surface depends on the local sea 
state, wind speed, and atmospheric stability. This study 
assumes that the drag coefficient is a function of height 
above the sea surface only. With this assumpt10nm yee 
variation of wind speed with height directly above the sea 
surface is a log profile, consistent with the "consi 
flux" approximation for the atmospheric mixed layer Jane 
Wind speed meaSurements for both data sets are for a height 
of 19.5m above the sea surface, where a constant drag 


coefficient of 0.0013 is assumed for use in the ocean model. 
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Garwood's (1977) ocean mixed layer model uses wind speeds at 
any height above the surface, but the drag coefficient needs 
to be adjusted approximately. Therefore, based on a log 
profile for the change in wind speed with height above the 
Surface the equivalent drag coefficient at 10 m is 0.0015. 
Firegmpiiets a ScCatterplot on a log=log scale of the 
Seuriace stresses a owilabee from the SASS and FNOC wind 
[Mmeelds, It is important to compare these fields directly 
because of the non-linear relationship between the wind 
speed and the surface stress. Eee 1 | Showsme a good 
Swereclatvion between UmewSuriace stress fields. The) scatter 
is relatively small for large values of the stress. fhe 
points are close to and relatively symmetrically distributed 
eels the line of perfect correlation. This indicates good 
overall agreement between the SASS and FNOC wind stress 
values, especially for average to high wind stress values. 
The scatter becomes larger for wind stress values below 0.01 
dynes/om?. However, this constitutes less than 1% of all 
the data points. Thus, these wind stress values derived 


from the SASS and FNOC fields are quite similar. 
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Figure 11. Seatterplot on log-log scale of surface  Sstigege 
Values computed from SASS and FNOC surface wind 
speed fields. The diagonal is included for 


bias detection. 
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lee ODMR COMPARISON 


A. PODEG@r ruorsics 

The model used in this study, described in detail by 
memewood (1977), 1s a second order closure, bulk turbulent 
mugeceme enersgy @IKE) model which separately retains the 
horizontal and vertical TKE equations. The model forcing 
feameloauosS Of the Surface fluxes of momentum and buoyancy. 
The surface momentum flux increases the horizontal TKE 
through the interaction of the surface stress with the 
vertical shear of the mean current. The buoyancy flux is 
produced by latent and sensible heat exchange, long wave 
Mmuddation at the Seles, Deavitnemmeys SOlar radiation 
throughout the euphotic zone and the exchange of heat at the 
mixed layer's lower boundary by entrainment. The model 
assumes that the vertical turbulent velocity (sealed from 
the vertical TKE) within the mixed layer transports the 
Surface generated TKE to the entrainment zone. The rate of 
entrainment of cooler water within and below the thermocline 
1S proportional to the amount of vertical TKE within the 
ayer, Pie etn edie oreo tinea varie redistribution 
mor izontal TKE from pressure rate of strain action, or by 
surface cooling, and it is decreased by dissipation and 


entrainment. The TKE available in the entrainment zone 
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erodes the underlying stable water mass, deepening the mixed 
layer. The mixed layer shallows when there is insufficient 
TKE to deepen or maintain the current mixed layecre@eaaa 
against the existing thermocline and a positive buoyancy 
Flux into Chewlaven 

Since the model depends on the interaction between the 
heat and momentum fluxes to determine the amount of TKE 
available for changing the mixed layer depth (MLD), it is 
important that the phase and magnitude of the surface wind 
forcing and surface heating be properly represented in the 
data sets used for model boundary conditions. AlSO =e ae. 
important that the data set used to initialize the model's 
thermal structure portray accurately the thermocline region 
below the mixed layer. The steepness of the temperature 
gradient in the thermocline region will determine the amoune 
of available TKE necessary to deepen the mixed layer a given 


aun Wier 


B. 30 DAY SIMULATION 

A 30 day simulation of the mixed layer was performed 
both with the SASS wind fields and with the FNOCG (uaa 
fields. The initial sea surface temperature fields were 
subtracted from the resulting 30-day Simulation fields to 
yield the temperature change fields shown in Fig. 12. 

The model predicted cooling over most of the region 


when the SASS winds were used and over about half the region 
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Figure 12. (a) 30 day model simulation of temperature change 
Gi) July 15 to August 15 at sea level using SASS 
WinGee Stress “boundary sconditiion: CeewO 5s 
(b) Same as (a) for FNOC wind fields. 
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when FNOC fields were uséd. This prediction of nev coomaae 
is unlikely for the July to August time frame in Ghee 
are of 1G. However, it is unlikely that the wind stress 
boundary condition was to blame, because both wind ia 
created similar trends. It is more probable |thavuaaie 
cooling was caused by some other forcing variable. Possibly 
suspect are the heat flux fields Since they are a 20g 
climatological average whereas the wind fields and the 
initial ocean thermal structure are not. AlSO SWS Qe 
the initialization field which may not have sSuffiecia, 
vertical resolution to adequately depict the deatea 
gradient of the thermocline region. This underestimate of 
Che initial tangential gradient would decrease the amountwon 
TKE needed to deepen the mixed layer, which would Cause 
MLD to inerease too much, lowering the mixed Sie 
temperature. A preliminary investigation indicates that the 
initial temperature structure fields may be the} praia 
cause of the unseasonable cooling simulated. 

Regardless of the reason for the cooling trend aaa 
wind fields may still be compared by way of model output 
because all other model variables remained unchanged during 
the simulations. Fig. i2 shows that the wind fields progime 
similar patterns. The maximum cooling and warming takes 
place in the same regions. However, these resultsS Show a 


sensitivity to changes in the Wind Stress bound 


2 


mir olhotionular Was Unexpected in light of the seemingly 


small differences in Wind speed shown in Fig. /f. 


on TIME oohn thos Or SO VAY SIMUBATION 

To further investigate the reasons for the differences 
in the 30 day temperature change fields between SASS and 
Pee forcing, three locations were chosen for time series 
comparisons. 

1. Location 36N, 180W 

This location was chosen to coincide with a warming 
ma eeron in born Simulated Pields. Fig. 13 18 a plot of the 
Mme series Of Wind Speed and daily maximum mixed layer 
@apcon, for both SASS and FNOC forcing. 

PrOtm@eune Tresulesepresentedein Fig. 13, the model's 
sensitivity to wind speed differences as large as 4 m/s is 
fomeia. Fig. 13 illustrates that once the MLD predictions 
Meeerer Significantly due-to differences in the wind Speed, 
any subsequent agreement between the MLD forecasts is likely 
to be only coincidental. This results because the forecast 
emeene MLD change 1S dependent on the current MLD. As the 
MLD grows, it erodes the thermal structure below the mixed 
[Meer , Which, in turn, reduces the rate of further 
deepening. THLlS iS beste lUsSerarvedsduring the five days 
moelLOwWwing, Day 17. During this time, the wind speeds 
merrelate very well, but the changes in MLD do not. The 


meeps Drceceeding Day if SUifieclently altered the thermal 
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(a) Time series of SASS (solid) and FNOC (dotted) 
Wind speeds during 30 day simulation at 36N, 180. 
(b) Same as (a) for maximum mixed layer depth 
per 24 hour period. 
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acto SO tMal the subsequently similar Wind stress 
foreing produced very different changes in MLD. 
2. Location HON, 175W 

This location was chosen because it was in the 
region where the SASS and FNOC isotach fields shown in the 
Mmoependix differed the most over the 30 day period. The 
results in Fig. 14 show model sensitivities similar to those 
mound in Fig. 13. ite Siwrsteiour days siow a verys good 
agreement between the wind fields and the resulting phase 
agreement between the simulated MLD's. Although the wind 
swoedso diverged sharply for short time periods, sometimes 
cauSing large differences between the fields, the MLD 
Simulations remained remarkably well in phase throughout tne 
mumme series. This result indicates that the duration of the 
difference in the wind speeds is an important consideration 
in the model's sensitivity to these changes. 

3. Location 42N, 155W 

lms whOCAbLL On NdasmmcnOosSem toe COlncide witt the 
fmeero0n Of Maximum cooling on both forecast fields. The 
Moree. ts at this location are shown in Fig. 15. ies 
Gemolbts also indicate the importance of a consistently 
accurate wind speed measurement for ocean mixed layer 
meacling. Once the wind speed measurements differ, the 
thermal structure in and below the mixed layer is altered 


and no subsequent correlation between the wind speeds is 
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likely to bring the simulated MLD changes into magnitude and 


phase agreement. 


De DISCUSSIONSOR PR ESU rs 

Primarily, these results show that the model @seuew 
sensitive to the variations in the wind speed boundary 
condition. Although the wind fields were Similar by direc? 
comparison, the differences wnen used as model forcing are 
Significant. This illustrates the importance of a@euiaaee 
Wind speed measurements for ocean mixed layer model use. 
The results also indicate the importance of NOU Samme ae 
Wind speed accuracy, but also the consistency (Gia 
accuracy throughout the forecast period. The model will 
maintain a "record" of these wind speed inaccuracicesmauaama 
e€ach subsequent Cime step of tne simubacvoneat tera have 
been introduced because of their effect on the thermal 
Simeulue Qure. Thus, differences appear to be cumulative, 
causing divergence in the mixed layer calculations. 

The results demonstrate that the wind speed boundary 
condition in this study is the primary forcing me@@iiaeaee 
influencing tne Synoptie scale Changes nee Thales 
result, however, should be expected since the climatological 
heat flux boundary conditions used in this study do nov 
contain Synoptic scale Gite rua vi ens. Additionally, @yae 
inaccuracies appeared to be less influential on MLD changes 


when they altered deepening rates than when they altered 
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shallowing rates. This may be due to the model's technique 
for deepening, which is done in several smaller time 
macremenvs according to the TKE available during the time 
step for mixing, aS compared to that used for shallowing, 
which is done in one calculation for the time step. 
[MeremorewremMme GensiStent accuracy Of Synoptie scale 
Winds used for forcing ocean mixed layer models is an 
extremely important factor in accurately predicting synoptic 
feaere Ghanges in MLD. ~The results of this study imply the 
importance of accurate surface heat fluxes. See Jaramillo 
(1984) for a more complete discussion. ies ene ae 
dependence of ocean synoptic seale MLD predictions on both 
Mmremewind Stress boundary condition and the current MLD 
meee aces ene importance of the accuracy and resolution of 


Moe Initial thermal structure fields. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY OF SRE Suis 

It has been shown that with quality (coi 
restrictions carefully imposed on the SASS data SeGuyyyea- 
SASS wind measurements compare very well in pnase and 
Magnitude with the FNOC wind measurements over a large ocean 
region on a grid which resolves tne synoptic scale vue. 
primary cause of poor comparisons were the occasionally 
large data gaps in the SASS wind fields which were filled by 
IOAN aoe WENy, Odie These gaps were largely attributable to 
changes in the SASS operational mode during, 
acquisSivuilome. A secondary cause of weak comparison was the 
inability of the FNOC wind fields to resolve events on 
Spatial scales as small as those resolved by the SASS wand 
fields. There was shown to be no significant differen@oman 
bias between the wind stress values derived from the SASS 
and FNOC fields. A comparison of ocean model simulations 
using the two wind fields in Garwood's model (1977) revealed 
a high sensitivity by the model to any Variatlonsmae 


SyYNOptie scale Wind Varese 
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Be CONCLUSIONS 

Iitsesemeeweias GComparedstne SASS and FNOC wind fields 
ea uci r Umit y Tor oceamemixed dayersemodeling. These 
fields have very different possible error sources, stemming 
from their respective methods of data acquisition and 
pealysis. Even™after imposing quality control restrictions, 
the SASS wind fields suffered from several factors: the use 
Soe an Cmpirical algorithm to derive wind speed, the 
mmaeeecuracy Of interpolating to fill data gaps, and 
differences between the magnitudes of aliased wind vector 
measurements. The FNOC wind field accuracy suffers from a 
lack of sufficient numbers of observations over large ocean 
regions, compounded by the limitations of atmospheric model 
Meores introduced by simplifying assumptions. In light of 
mimese Widely differing possible sources of error, it is 
Significant that these fields agree as well as they do, 
mmoamng credibility to both methods. In particular, the 
verification by the FNOC wind fields of the accuracy of the 
SASS wind measurements when averaged spatially and 
temporally indicates tne potential of the secatterometer wind 
measurements for ready use in any ocean mixed layer model 
which does not require wind direction. 

The ocean mixed layer model used in this study was very 
sensitive to the variations of the wind fields used as 
boundary conditions. The accuracy required by the model is 


higher than the verifiable accuracy of either wind field. 
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However, though the sensitivity of the model to the wind 
field has been shown, it should be noted that monthly 
forecast accuracy is limited even more by difficulties@in 
specifying the heat flux boundary condition and )tnaaaa- 
initialization fields. The SASS fields do providewiaea 
fields which are comparable to the FNOC wind fields and 
which may be used as an alternative data source for further 
study of the ocean mixed layer model's response to wind 
here ae. 

Perhaps the most promising aspect of Uhis (sua 
results is that it creates a basis for combined use "Ogayew 
and FNOC wind fields. While studying the im pagieuow 
scatterometer data on weather prediction model forecasts, 
Atlas et al. (1982) referred to this potential im @ieaeg 
COme lus wonicn 

"Regional studies demonstrate that Synoptic 
subsynoptic systems which are not detected Or “pG@@mam 
analyzed are accurately represented in scattCerometer data, 
and short-range forecasts based on this data should be 
improved," 

While SASS has the potential to significantly improve 
the detection and analysis, of subsynoptic and Syne@puae 
systems, the use of a dynamic model to fill gaps in thevsre 
coverage would be a substantial improvement over non- 
selective use of polynomial interpolation. Whe potencam 
for combining the SASS data and standard atmospienmme 


analyses 1s promising for ocean mixed layer model 


applications because of the symmetry of Che Winds 
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meaeverplove comparison. indicating that both wind fields 
Pmogdguce COmparable wind stress values. 

Additionally, the scatterometer scheduled to be on 
board the NROSS satellite should increase the potential of 
Mes data source (Satellite Surface Stress Working Group, 
mge2). Tavis OCmu=eaeapre Or reducing the directional 
ambiguity in measured wind vectors, as well as providing a 
Perel eseolurion Of 25kmom— ALSO, should 1t complete its 
three year operational mission, it will provide improved 
data coverage by operating in a double side swath mode 
continuously and by increasing and stabilizing the average 
heeaeh Of Side Swaths. Operation in the double side Swath 
mode could easily reduce the time period required to produce 
Pee ate coverage for interpolation to l2 nours or less. 
Furthermore, FNOC now uses a more advanced operational 
atmospheric model, the Naval Operational Global Atmospheric 
Peediction System (NOGAPS), which may improve both the 
coverage and quality of both the surface winds and the 


surface heat fluxes over the oceans. 


Oe RECOMMENDATIONS 

Hoe iscmrccommended thav further study using the SASS 
data set, available from JPL, with Garwood's ocean mixed 
layer model (1977) concentrate efforts on three areas. The 
minst 1S am extension to this study, which explores in more 


aoe Oomueodch S SCUSIUlLViLy COyUNe Synoptic variations of 
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the wind speed boundary condition. The SASS fields have the 
potential for Synoptic mixed layer modeling use in related 
areas of the reduced ADS region. A more formal SUatGaisae 
analysis of the SASS wind fields using time-space objective 
analysis methods may enhance this potential. The second 
effort is to develop an acceptable method for combining the 
SASS and FNOC fields and using the combined field for model 
forcing. Finally, using SASS and FNOC wind fields, age 
of the relationship between the model's response to synoptic 
and mesoscale forcing 1S in order. An investigation (ota 
potential of using local synoptic forecasts combinegaamma 
regional mesoscale forecasts may aid sine this seu 

Should the NROSS scatterometer make longer time period 
data sets available, a similar comparison study using a 
large number of monthly forecasts should be mMademera 
determine quantitatively the differences in the model's 
response to the scatterometer and FNOC (NOGAPS) wind fields, 
including seasonal and regional changes. The NROSS wind 
fields should allow comparison with a l2 hour tineeeeoeeee 
such studies will require not only NROSS and NOGAPS 7 tea 
for model forcing, but also ocean thermal structure data for 
model initialization and ver ist7 ca tumeme This and Oia 
Studies (e.g. Jaramillo, 1984; Stringer, 1984) have shown 
Chat the ADS NORPAX analysis is barely adequate. Tits 


future large scale ocean thermal structure analyses will 
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have to be improved if ocean model forecasts are to be 


improved. 
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APPENDIX A 


SASS AND FNOC COMPARISON 


The following pages present figures which compare SASS 
data coverage before interpolation and SASS and FNOC isotach 
analyses for the 34 day period of this study. The figures 
are Similar to Fig. 11 and are arranged four to a page 
Table A-1 lists the date corresponding with day ShOwnwiaen 


Cach ti eure. 
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List of Time Periods Studied 


JULIAN DAY 
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